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Abstract 

The role of the Taub time gauge in cosmology is linked to the use of 
the densitized lapse function instead of the lapse function in the varia- 
tional principle approach to the Einstein equations. The spatial metric 
variational equations then become the Ricci evolution equations, which 
are then supplemented by the Einstein constraints which result from the 
variation with respect to the densitized lapse and the usual shift vector 
field. In those spatially homogeneous cases where the least disconnect 
occurs between the general theory and the restricted symmetry scenario, 
the recent adjustment of the conformal approach to solving the initial 
value problem resulting from densitized lapse considerations is seen to be 
inherent in the use of symmetry-adapted metric variables. The minimal 
distortion shift vector field is a natural vector potential for the new York 
thin sandwich initial data approach to the constraints, which in this case 
corresponds to the diagonal spatial metric gauge. For generic spacetimes, 
the new approach suggests defining a new minimal distortion shift gauge 
which agrees with the old gauge in the Taub time gauge, but which also 
makes its defining differential equation agree with the vector potential 
equation for solving the supermomentum constraint in any time gauge. 



1 Introduction 

In the 1950s, pioneering work by Lichnerowicz P^, Choquet-Bruhat [2] and 
Dirac |3] led to the Hamiltonian analysis of general relativity by Arnowitt, 
Deser and Misner (ADM) 0] using their ADM decomposition of the spacetime 
metric variables ( ^g a p) into a spatial metric (g a b) on a spacelike slicing of 
spacetime threaded by time coordinate lines whose tangent determines the lapse 
function (N) and shift vector field (/3 a ) variables, terminology introduced by 
Wheeler More recently the densitized lapse a = N/g 1 / 2 , namely the lapse 
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function divided by the square root of the spatial metric determinant (a spatial 
density of weight —1), has emerged to play a key role in better understanding 
the gravitational field dynamics and constraints in general relativity, especially 
in work by Choquet-Bruhat and York and their collaborators during the past 
decade (see 011] for a complete list of references, including Teitelboim [El El) 
Ashtekar JU| and Frittelli JH). In addition to the name densitized lapse, this 
function has been called the lapse density, slicing density, slicing function and 
lapse antidensity (because of the weight —1). York has suggested calling it 
simply the Taub function in honor of the man who put it to its first significant 
use in general relativity at the very beginning of the 1950s ^2]. We adopt this 
name here and follow York's new convention of letting a denote it instead of 
the lapse function TV as in his older work. 

These new variables of the 3 + 1 decomposition of the spacetime metric repre- 
sent local changes of the metric variables which lead to new linear combinations 
of the field equations in the Lagrangian or Hamiltonian approach. The original 
spacetime metric variables ^g a /3 lead to the field equations in (contravariant) 
Einstein form, equating the Einstein tensor to the gravitational constant times 
the matter energy-momentum tensor 

(4) G a/3 _ kT «/3 _ ^ 

Changing to the lapse, shift and spatial metric variables in the variational ap- 
proach automatically splits these field equations orthogonally into the spatially 
projected Einstein equations (the "evolution equations" ) and the Einstein con- 
straints. The further change of variables to the Taub function, shift and spatial 
metric variables reshuffles these equations again to yield the evolution equations 
in Ricci form, namely the spatial projection of the Ricci form of the Einstein 
equations 

Wr«0 = K (T afs - - ^g ap T\) , (2) 

and the same Einstein constraints. 

This mixed form of the Einstein equations (Ricci evolution equations plus 
Einstein constraints) has long been known to be rather useful for studying var- 
ious problems in general relativity. It was used by Darmois |13| during the 
period from the 1920s to the 1930s and continued by his student Lichnerowicz 
starting in the 1940s and by Choquet-Bruhat in the 1950s and beyond, in- 
cluding York in his review article on the kinematics and dynamics of general 
relativity in 1979 Most recently it has been recognized that with the Ricci 
evolution equations combined together with the additional evolution equations 
for the constraint functions coming from the Bianchi identities and source field 
equations one has a symmetric hyperbolic system for the constraints with pos- 
sibly important consequences for numerical solution since the constraints then 
propagate in a causal domain of dependence in contrast with the corresponding 
system using the Einstein evolution equations |11|." 

The Taub function indirectly made its appearance in Taub's classic paper 
[TB*] introducing Bianchi cosmology in 1951, where he finds his Taub solution of 
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the diagonal vacuum Bianchi type II case in the time gauge N — g 1 ^ 2 , namely 
that the Taub function is unity: a = 1. Misner JHj used this same time gauge 
(see his Eq. (44), recalling that e~ n = g 1 ^) which he called the supertime time 
gauge, although a slight (affinely related) variation of it a = 12 is more appro- 
priate in the full Hamiltonian analysis since it makes the kinetic energy term 
in the gravitational Hamiltonian for the diagonal Bianchi models correspond 
exactly equal to the one associated with the 3-dimensional Minkowski met- 
ric in terms of the natural logarithmic isotropy/anisotropy variables (3° = — SI 
and which parametrize the diagonal metric coefficients. This generalizes 
nicely with the DeWitt inner product on the space of inner products (or simply 
DeWitt metric ^7j) to higher-dimensional spatially homogeneous spacetimes 
|18| . Looking back decades later, the Misner (and later Ryan ^|) Hamiltonian 
analysis is a bit hard to follow since they always start with the reduced Hamil- 
tonian in the O time gauge for which a — N/g 1 / 2 = 12(p ) _1 = 3/(g 1 / 2 TrX), 
where H = |po| is the (reduced) Hamiltonian, and po is the momentum dual 
to the isotropy variable f3° = —ft = ^ln<?, rather than starting with the La- 
grangian/Hamiltonian of the general theory and reducing the latter. In the case 
of the Taub Bianchi type II vacuuum solution, the Taub time gauge leads to 
the Hamiltonian problem of scattering in a time-dependent 1-dimensional ex- 
ponential potential in flat 3-dimensional Minkowski spacetime, which is easily 
decoupled by a Lorentz transformation into three independent 1-dimensional 
exponential scattering problems, one with vanishing potential corresponding to 
free motion, the others easily solved in terms of hyperbolic functions 20 21 22 . 
For Taub's locally rotationally symmetric Bianchi type IX vacuum solution, the 
free motion degree of freedom is suppressed by the local rotational symmetry. 

For spatially homogeneous spacetimes of any dimension, the term Taub time 
gauge has been well established to describe the lapse condition a — constant, 
i.e., the Taub function is fixed to be a constant. It seems reasonable to keep the 
terminology Taub time gauge to refer to a spacetime constant Taub function 
function and let the "generalized Taub time gauge" (or some suitable new name 
yet to be determined) describe the condition that a be time-independent in a 
general inhomogeneous spacetime, i.e., the Taub function is independent of the 
time coordinate in the associated coordinate system. Of course this includes 
the Taub time gauge as a special case. Choquet-Bruhat and Ruggeri found that 
exactly this latter generalized Taub time gauge converts the 3 + 1 split Einstein 
equations into a hyperbolic system j2Hl when expressed as a second order system 
in the extrinsic curvature. They noted that when the shift vector field is zero, 
the generalized Taub time gauge corresponds to a harmonic time coordinate 
slicing, namely t. a ' a — 0. A short calculation of the spacetime Laplacian of t in 
3 + 1 variables 

t a ' a = N-^dtXna + NVaiN- 1 ^)] 

= a^N-^dt - £p\a (3) 

where the Lie derivative of the weight —1 density a [23] is 



£pa = f3 a V a a - aV a /3' 



(4) 
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(which corrects (13.24) of EH| and corresponds to Eqs. (22)-(25) of HJ with 
/ = 0) shows that the vanishing of t- Q ' a is equivalent to the Taub function being 
time-independent when the shift vector field (3 a is zero (or more precisely, when 
the relative velocity N^ 1 ^ between the time coordinate hypersurface normal 
direction and the 4-velocity tangent to the time coordinate lines is divergence- 
free in the spatial geometry). This is exactly the case in the Bianchi type II 
and IX vacuum spacetime solutions of the Einstein equations found by Taub 
in zero-shift gauge, where the Taub function, already spatially constant as a 
spatially homogeneous scalar, is a spacetime constant in harmonic time gauge, 
the gauge he used to solve the equations. 

The Taub time gauge was also used by Belinski, Khalatnikov and Lifshitz 
|23 ( see their Eq. (3.4)) in their analysis of the change of Kasner epoch due 
to a single curvature wall, the Bianchi type II approximation to the spatially 
homogeneous dynamics in the BKL limit approaching the initial singularity 
in a generic spacetime in which the dynamics becomes pointwise like spatially 
homogeneous dynamics. It should be noted that they always considered the 
Ricci form of the spatial Einstein equations in their analyses [2B] and they 
and others in the Russian school at the time (see for example) often used 
the Taub time gauge to state Einstein's equations for Bianchi models in the 
form of the Ricci evolution equations plus the Einstein constraints. For the 
diagonal such models they considered, the time derivative terms in the Ricci 
evolution equations in the Taub time gauge reduce to second time derivatives of 
the logarithms of the diagonal spatial metric coefficients, the simplest form they 
can take, corresponding to the flattened DeWitt metric (having scaled away its 
spatial metric determinant factor) on the space of these logarithmic variables. 

The Taub time gauge considerably simplifies the evolution equations by elim- 
inating the extra terms involving the kinetic energy as a multiplier that arise 
from the differentiation of the factor g' 1 ^ 2 in the kinetic energy part of the 
Hamiltonian. A calculation equivalent to the explicit vacuum calculation of 
Anderson and York (SO) (see their Eq. (12)), shows that this is equivalent to ob- 
taining new evolution equations which are only linear combinations of the spatial 
projection of the Ricci form of the Einstein equations, even in the nonvacuum 
case. Jantzen [20112] (see also ^Hl) made this claim without showing the cal- 
culation, and used the Taub time gauge and Ricci form of the field equations to 
explain some exact Bianchi cosmological solutions which are easily obtained in 
the Taub time gauge, for which the evolution equations (see his Eq. (8) in [21]) 
are particularly simple for a diagonal metric. These two articles both promote 
the use of the mixed form of the Einstein equations: Ricci form of the evolu- 
tion equations plus the Einstein constraints, an approach which was exploited 
in |31| for interpreting vacuum solutions of the Ricci evolution equations for 
which certain multiples of the Einstein constraints are then conserved, leading 
to nonvacuum solutions of the full equations. 

The Taub function proved very useful in investigating the dynamics of the 
Bianchi cosmological models. The variable x = a -1 , namely the reciprocal of 
the Taub function, was introduced already in Jantzen |31j as a convenient vari- 
able which in the Silkos time gauge N — (533) 1 / 2 (equivalently a = (511522) 1 ^ 2 ) 
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for certain Bianchi types and initial data has simple exponential or hyperbolic 
function solutions. Later the same symbol x^ 1 was explicitly used for the Taub 
function due to the simplifying properties of the Taub time gauge and Ricci 
evolution equations in the spatially homogeneous dynamics |S2 El- In- 
terestingly enough Anderson, Choquet-Bruhat and York also initially used the 
reciprocal Taub function which they also called a [33] (see their Eq. (66)) before 
switching to the Taub function itself. Note that a was the symbol originally 
used by York for the lapse itself. 

This same simplifying feature of the Taub time gauge was used by Moncrief 
|36| in his study of Einstein equations for spacetimes with one Killing vector, 
although the gauge applies to the Weyl rescaled metric in his fiber bundle fomu- 
lation so it is not immediately apparent, as noted explicitly in |22| (see section 
VI). In fact the more general power law lapse conditions, especially in higher- 
dimensional spacetimes involving dimensional reductions, together with other 
variable transformations, seem to be the larger setting in which the simplifica- 
tions of the Taub time gauge extend to other special symmetries. 

Teitelboim introduced the Taub function in his study of the Hamiltonian 
constraints (see his Eqs. (4.1) and (4.4) in [H], where he uses the symbol N 
for the Taub function function) to simplify considerably the Feynmann path 
integral approach to quantization of the gravitational field, since it removes the 
square root of the spatial metric determinant from the DeWitt metric on the 
space of spatial metrics in the Hamiltonian formulation of general relativity, 
simplifying the path integral measure considerably. He calls the constant Taub 
function time gauge together with zero shift vector field the "proper time gauge" 
(see his Eq. (4.20a) in |S]), which is a clear misnomer since this is not associated 
with a proper time coordinate gauge, acknowledged just before his Eq. (2.2) in 

but arose from the analogy with the dynamics of a relativistic point particle. 
Apparently similar considerations motivated Ashtekar to introduce the Taub 
function in his "new variables" (see his footnote 17). 

Section 2 shows how two simple changes of variables, first from the space- 
time metric coefficient fields to the lapse, shift and spatial metric and then to 
the Taub function, shift and spatial metric associated with a given slicing and 
threading of spacetime, repackage the Einstein equations in terms of constraints 
and evolution equations when derived from a variational principle. Section 3 
shows how these changes in evolution equations affect the evolution of the con- 
straints. Section 4 reviews the final version of the conformal approach to the 
initial value problem of solving the supermomentum constraint using a vector 
potential, which has been improved by Pfieffer and York Z7 based on Taub func- 
tion considerations, and suggests extending this improvement to the minimal 
distortion shift gauge so that its defining differential equation agrees with the 
vector potential equation. Section 5 shows how the natural symmetry-adapted 
decomposition of the metric variables in spatially homogeneous cosmology nicely 
reflects the new approach. 
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2 Local variable changes and the chain rule in 
the variational approach 

Starting with a timelike foliation of the spacetime and adapted coordinates 
{x a } = {x° = t, x a }, the spacetime metric line element 

ds 2 = {i) g af 3dx a dx f3 = -N 2 dt 2 + g ab (dx a + f3 a dt)(dx b + (3 b dt) (5) 

can be re-expressed in terms of the ADM variables: the lapse function N, the 
shift vector field (3 a and the spatial metric g ab (inverse denoted by g ab ). This 
can be thought of as a local change of metric variables 

( (4 W (4) 3o a , (4 W) - (-(N 2 - 3afc/3 a /3"), g ab p\ 9ab) -► (N, f3 a ,g ab ) , (6) 

where the absolute value of the spacetime and spatial metric determinants ^g — 
I det(g a p)\ and g = | det(g a b)\ are related by the square of the lapse ^g — N 2 g. 
The associated differential is 

(d^g 0Q ,d(% 0a ,d(% ab ) 

= (-2NdN + 2(3 a dp a + P a p b dg ab , g ab df3 b + S^ c a S d \(3 b dg cd , dg ab ) , (7) 

while the contravariant metric variables are 

( ( V°, (4) 5° a , {A V b ) = {-N~ 2 , N~ 2 (3 a ,g ab - N- 2 B a 3 b ) . (8) 

With the standard variational formulas 

*<V /a = iWf^, 
f<V = - ( V 7(4) /^ (4 W , 

S{ (4) 5 V2 (4) gQ /3 (4)^) = (4) ff l/2 ( 1 Wg cry (A) g f)6 (4)^ _ (4)^^ (4)^ 

+ (4 V /2( V^ (4) i? a/3 , (9) 

the usual scalar curvature gravitational Lagrangian produces the lefthand side 
of the field equations in Einstein tensor form upon variation 

L = <V /2 (4) ^ ,6L = - {i W 2 (i) G a Pdg af} + BTs , (10) 

where "BTs" (boundary terms) stands for additional divergence terms which in- 
tegrate out to the boundary of a suitable tubular integration region between two 
time coordinate slices compatible with the slicing and its adapted coordinates, 
where variations are fixed and so make no contribution to the field equations 
under suitable conditions on the metric variables. Of course these boundary 
terms are themselves important for various reasons, but not for the present dis- 
cussion and so their details will be omitted. York pointed them out in the 3 + 1 
approach in 1972 |2Hj, while soon after they were analyzed from a different point 
of view by Regge and Teitelboim jSH]. Wald has discussed them in appendix E 
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of his textbook |JDj . One must add terms to the various variational principles 
to cancel out unwanted boundary terms on the bounding hypersurface of the 
spacetime region of variation, since only the metric variables themselves can be 
held fixed there and not their off-hypersurface derivatives. 

This Hilbert Lagrangian in turn differs from the ADM gravitational La- 
grangian only by such boundary terms Ladm = L + BTs and so produces the 
same left hand side of the field equations 

5L ADM = -Ng^G ab S ^g ab + 2 ^G ob S ^g Qb + ^G°°5 ( 4 > 5oo ] + BTs . (11) 

Explicitly 

Ladm = Ng 1 ^ 2 [K a c K c b — K a a K b b + R] 

= n ab g ab — NTi — f3 a H a + BTs , (12) 

where / = df/dt is the ordinary time derivative acting on a function (equiv- 
alent to £o t when acting on tensors since the coordinates are dragged along), 
and the extrinsic curvature, canonically conjugate momentum, supermomentum 
and super-Hamiltonian are, letting a vertical bar denote the spatial covariant 
derivative X a \ b = V b X a ), 

K ab = -^[9ab - £ fi 9ab] , 7T ab = dL ADM /dg ab = -g^ 2 [K ab - g ab K c c ] , 

H a = 2Ng^ ^G° a = 2g^ 2 = -2g^(K b a - K c c 5 b a )\ b = -2* h a \ b , 

U = -27VV /2(4) G 00 = 2.g 1 / 2 ( 4 )G ± ± 

= g l ' 2 [K a c K\ - K\K b b -R] = g' 1 ' 2 ^ C n c b - ^ q tt\ - R] . (13) 

Here the "perp" indices refer to components along the timelike unit normal 
vector field to the spacelike slicing e± = N^ 1 (dt — (3 a d a ) and its sign-reversed 
index-lowered 1-form 8 1 " = Ndt; note that K ab = — ^£ e± g a b- The first two 
terms in the super-Hamiltonian are the kinetic energy terms corresponding to 
the DeWitt metric [T7] . 

Simply re-expressing the variation of the ADM Lagrangian using the chain 
rule for the change to ADM metric variables automatically leads to new lin- 
ear combinations of the field equations which reflect the spacetime splitting of 
the symmetric second rank tensor into a space-space projection, a time-space 
projection and a time-time projection under the orthogonal decomposition of 
the tangent space into a spatial subspace tangent to the slicing and its time- 
like normal direction. Contravariant indexed components correspond to the 
time projection (since = is the sign-reversed index-lowered normal), 

while covariant spatial indices correspond to the spatial projection (since e a are 
tangent to the spatial hypersurface of constant time) 

SL ADM = -Ng 1 / 2 [g ac g bd ^G cd Sg ab + 2^G° b Sp b ~2N^G m SN + BTs] 
= [-Ng 1 l 2 g ac g bd ^G cd ]Sg ab + H b 5(3 b + H5N + BTs . (14) 
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The spatial projection of the Einstein tensor leads to the left hand side of the 
evolution equations in Einstein (tensor) form, while the remaining field equation 
left hand sides are the gravitational super-momentum and super-Hamiltonian 
constraints, namely the Einstein gravitational constraint functions. (Obviously 
the time-space projection of either the Einstein tensor or Ricci tensor form of 
the field equations produces the same supermomentum constraint.) 

A final change of variables expressing the lapse function in terms of the Taub 
function a = N/g 1 / 2 , with differential 

dN = diag 1 ' 2 ) = g x / 2 da + ag 1 ' 2 g ab dg ab , (15) 

leads to a final shuffling of the field equations 

6L ADM = -Ng 1 / 2 [g ac g bd ^G cd + g ab ^G ± ± }Sg ab 

+n b S/3 b + g 1 ' 2 U5a + BTs . (16) 

The new left hand sides of the gravitational evolution equations are now in a 
trace-modified Ricci tensor form, as a side calculation shows. 
The orthogonal splitting of the trace of the Ricci tensor is 

W R = _7V2 (4 )i? 00 + g cd (4) Rcd = (4) R ± ± + g cd (4) Rcd ^ 

which implies that 

(4) G ± ± = _ 7v 2(4) G 00^ ( _ Ar 2 )( (4) i? 00_ l(4) ff 00(4 )i?) 

= l( {i) R ± ±-9 cd{i) R cd ) (18) 

and so 

Wr ± ± = 2 ^G ± ± + g cd ^R cd . (19) 
Then the spatial projection of the Einstein tensor satisfies 

(4) G ab + g ab ^G\ = ^R ab \g ab ( «W ± + g cd ^R cd 2 ^G\) 
= (i) Ra b -9a b g cd{A) R cd 

= (S c a S d b -g ab g cd )^R cd = S ab . (20) 

This transformation associated with index-shifting by the DeWitt metric on the 
space of symmetric tensors only changes the trace of the spatial projection of 
the Ricci tensor by the factor 1 — 3 = — 2 and is easily inverted by 

W Rab = (8 c J d b - \g ab g cd )S cd , (21) 
leading to the final formula 



5L ADM = -Ng 1 l 2 [g ac g bd ~g ah g cd ]^R cd 5g ab 

+n b S/3 b + g 1 ' 2 U5oi + BTs . (22) 
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Thus with the new independent variables, the evolution equations are equivalent 
to the Ricci tensor form of those evolution equations, as shown by Frittelli JT] 
and by Anderson and York |41|. 

In spatially homogeneous cosmology, the Taub time gauge makes a constant, 
so imposing this time gauge on the Lagrangian automatically yields the Ricci 
evolution equations as the variational equations for the spatial metric. The sim- 
plification in the variation which produces these equations is obvious because the 
g 1 / 2 factor in the kinetic energy is removed, and no longer generates additional 
terms proportional to the kinetic energy (and hence to the super-Hamiltonian 
constraint function when combined with the corresponding term arising from 
the potential energy term, namely the spatial scalar curvature term). 

3 Evolution of the Constraints 

Working with the original coordinate frame shows clearly the simple changes of 
variable and their subsequent effect on the form of the variational equations. 
However, the orthogonal decomposition of tensors is very cumbersome since the 
frame is not adapted to the spacetime splitting parallel and perpendicular to 
the time coordinate hypersurfaces. By replacing the time coordinate derivative 
in the coordinate frame by its normal component, one obtains a partially or- 
thogonal frame adapted to this splitting which makes it trivial to split tensor 
equations 

e =d t - N°d a = Ne ± , e a = d a , 6° = dt = N^d 1 - , d a = dx a + N a dt . (23) 

This is called the Cauchy adapted frame for the foliation by Choquet-Bruhat 
0. Normalizing the normal and its dual one- form is even more convenient for 
calculation. 

One can generalize these frames to adapted computational frames |14| by 
letting e a be any spatial frame (each frame vector tangent to the time slices) 
which is Lie dragged along the time coordinate lines 

[e„,e 6 ] = C c ab e c , [d t ,e a }=0. (24) 

This is necessary in order to apply the present discussion to symmetry adapted 
frames in the case of spatially homogeneous spacetimes, for example. 

The behavior of the constraints depends on which evolution equations one 
uses. One splits the twice contracted Bianchi identities (divergence of Einstein 
equals zero) by evaluating components in one of the above adapted frames. Let 

M af3 = {A)Qaf3 _ Xhen 

M a ± , a = g- 1 / 2 d^g l ' 2 M^ 1 _) + N- l {NM c 1 _) ]c + M c d K d c , 

M a a , a = g- 1 / 2 £ ex (g 1 / 2 M x a )-N- 1 N {a M ± ± + N- 1 {NM c a ) ]c . (25) 

Provided that the matter equations of motion are satisfied, the divergence of 
the energy-momentum tensor is zero and hence these expressions must vanish. 
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In terms of the total super-Hamiltonian and supermomentum 

w (total) = 2 5 1 /2 M -L i i W (total) = 2g 1 / 2 M ± a , (26) 

these become 

g- 1/2 d ± (H (total) ) = -N- 1 {Ng~ 1/2 rC {tot&X)c )\ c -2M c d K d c , 
g- 1/2 £ e± (H a total) ) = N- 1 g- 1 ^N ]a H ( - taM) ~ 2N-\NM c a ) lc . (27) 
When the Einstein evolution equations M ab = hold, explicitly 

Ng V2 M ab = ji^ab + 2Ng -l/2^a c7r cb _}_ n ab n C c - ) 

-^^VWd ~ 7r c c 7r d d ) + Ng 1 / 2 {G ah - nT ab ) 

+Ng l ' 2 (N\ ab - g ab N c \ c ) = , (28) 

the constraint evolution equations become very simple and confirm that if the 
constraint functions are zero on the initial hypersurfacc, they continue to remain 
so during the evolution. 

To see the consequences for the constraints when instead imposing the Ricci 
evolution equations 

^R ab = n{T ab - ^g ab T) <- ^R ab - g ab ^R c c = K (T ab + g ab T^ ±) 

^0 = M ab + g ab M ± ± 

= £ eo n ab + 2Ng- 1 / 2 (iT a c TT cb - Vv,) 

+Ng l ' 2 {R ab - g ab R c c - nT ab - Kg ab T JL ± ) 

+Ng 1 / 2 {N\ ab - g ab N c \ c ) , (29) 

one must re-express the above equations in terms of the combination M ab + 
g a bM ± ±. Using the identity 

K = K C C = -^g ab £ e± g ab = -d^g 1 ' 2 , (30) 

one finds 

M a ±ia = g - 1 ' 2 d ± (g 1 / 2 M ± ± ) + N- 1 (NM c ± ) lc + M ± ± g- 1 / 2 d ± (g 1 / 2 ) 

+ [(M c d + S c d M^ ± )K d c ] , 
M a a;a = g- 1 / 2 £ e± (9 1/2 M ± a )-N- 1 N ]a M ± ± +N- 1 (-NM ± ± ) la 

+ [N- 1 (NM c a + NS c a M ± ± ) lc ] , (31) 

or 

M a Ua = ff- 1 a ± ( 5 M ± ± ) + JV- 1 (7VM c ± )| c 

+ [(M c d + 5 c d M ± ± )K d c ] , 
M a a;a = g- 1 / 2 £ e± (g 1 / 2 M ± a ) + N- 2 (-N 2 M ± ± ) la 

+ [N- 1 (NM c a + N5 c a M^ ± ) lc } . (32) 
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Again provided that the matter equations of motion are satisfied, one then 
has the result of Anderson and York 01] trivially extended to the nonvacuum 

case 

-2[(M c d + S c d M ± ± )K d e ] , 

g-^ 2 £ e± (nt tal) ) = Ar 2 (aV /2 W (total) )|« 

-2[N-\NM c a + NS c a M ± x ) lc ] . (33) 

Now when the Ricci evolution equations hold, the square bracketed terms are 
zero. However, as noted in the introduction, the combined Ricci evolution 
plus constraint equations supplemented by the propagation equations for the 
constraints themselves form a symmetric hyperbolic system, with the constraints 
propagating in a stable manner that is important for numerical solution, in 
contrast with the corresponding Einstein evolution system |11| . 

While the vacuum Einstein evolution equations are changed by the intro- 
duction of an energy-momentum tensor whose spatial part is nonzero, the Ricci 
evolution equations can remain unchanged if the spatial energy-momentum ten- 
sor is related to the energy-density by 

T ab + g ab T x ± =0 . (34) 

For a perfect fluid source, the energy-momentum tensor is 

T a/3 = pu a U(3 + p(u a up + (4) 5 Q/ 3) . (35) 

If this fluid is not tilted, i.e., is flowing orthogonally to the time hypersurfaces 
(4-velocity u a = n a — > u a — 0,u = 1), then T a f, = pg a b, T ± ± = —p, so if 
the fluid is stiff (p = p), then this condition is satisfied and the Ricci evolution 
equations are not changed. 

Thus if one has a solution of the vacuum Ricci evolution equations and 
the vacuum supermomentum constraint, one can use the super-Hamiltonian 
constraint function to define the energy-density p of a nontilted stiff perfect 
fluid when the gravitational super-Hamiltonian is negative, and it will satisfy 
d±(pg) — 0, leading to a conserved quantity for the equivalent nonvacuum 
system. 

This is the stuffing operation for spatially homogeneous and spatially self- 
similar spacetimes |31| . For such spacetimes one can also extend this re- inter- 
pretation of the vacuum Ricci evolution equation solutions with nonzero con- 
straint functions to the tilted case where the tilt is along an invariant direction, 
say e%. Decoupling and then solving the evolution equations requires that a 
be independent of 533, a convenient choice being a = (det g^) 1 ^ 2 , where g^ 
is the determinant of the remaining block of the spatial metric. One can then 
interpret solutions of the Ricci evolution equations with nonzero TL and H.3 as a 
tilted stiff perfect fluid. In the variational point of view, this time gauge leads 
to another set of evolution equations that are better suited to the particular 
symmetries of this class of models. 
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More generally the power law lapse time gauges [22] extend this kind of 
simplication to other scenarios in spatially homogeneous (and spatially self- 
similar) cosmology. In fact, almost all exact solutions of this type can be found 
by choosing the lapse function to be a power law function of appropriate metric 
components in order to obtain new evolution equations through the variational 
(or extended variational) approach which decouple for the appropriate choice of 
metric variables and thus admit solution |34j . 

4 The New Approach to the Initial Value Prob- 
lem 

To discuss the conformal approach to the initial value problem, we follow the 
notation of Pfeiffer and York [HZ]. Let an overbar denote the physical metric 
variables, which are related by a conformal rescaling to the unphysical variables 
without an overbar. The Taub function and the shift vector field and the trace 
t of the extrinsic curvature are not transformed 

a = a , P a = f) a , g ab = ^g ab , g 1 ' 2 = <fg x l 2 , N = 6 N , 
K\ = -^(g ac dg cb /8t - g ac £ g cb ) = -^^(g ac dg cb /dt - g ac £pg cb ) 

= A\ + ~r5 a b , (36) 

where A a b — K^ TF ^ a b is the tracefree part of the extrinsic curvature. The fixing 
of the Taub function under the conformal transformation is motivated by the 
nice properties of the Einstein equations that occur when this quantity is fixed. 

The mixed form of the Lie derivative term appearing in the extrinsic curva- 
ture has the transformation law 

9 ac £f39c b = g ac £pg cb + 41n^ c /3 c 5° 6 , (37) 
so its tracefree part is invariant 

[Lf3}\ = [g ac £pg cb } {TF) = [g ac £ P g cb ] (TF) = \ L PTb ■ (38) 

The same is true of the tracefree part of the time derivative term 

[g ac dg cb /dt]^ = [g ac dg cb /dt]^ , (39) 

suggesting that the conformal transformation of the tracefree part of the ex- 
trinsic curvature should be due entirely to the factor of the lapse. This is in 
fact reinforced by the conformal transformation properties of the divergence 
operator appearing in the supermomentum constraint. This constraint 



^b(K b a - K C c 8 b a ) = J, 



(40) 
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can be thought of as a condition on the divergence of the tracefree part of the 
extrinsic curvature, determining it in terms of the trace r and the source current 
1 =T X 

Ja 1 a 

V b A b a = ^V a r + j a . (41) 

To solve the supermomentum constraint in the Hamiltonian picture, namely 
in terms of the initial variables g ab , K ab and j a , one can conveniently choose a 
conformal representation of the tracefree part of the extrinsic curvature. Letting 
V a and V a be the spatial covariant derivatives associated with the two spatial 
metrics, the transformation of the divergence of a symmetric tensor S ab is easily 
evaluated 

gab = ^x-Agab Qr S a b = (j) x s a b 

-> VbS b a = ^ x (VbS b a + {x + Q)S b a Vb\n^-2S h bVa\n4>) . (42) 

Thus picking x = — 6 makes the divergence of a tracefree such tensor also 
transform by a conformal factor, namely by the weight —6 for the divergence 
1-form. This corresponds exactly to the transformation due to the reciprocal 
lapse factor above. One can decompose the tracefree part of any symmetric 
tensor into a transverse traceless part (zero divergence) and a pure divergence 
part using the tracefree Lie derivative operator. However, the covariant form of 
the divergence of the tracefree Lie derivative operator is invariant (if the vector 
doing the derivative is invariant) , so one must include an additional transforming 
factor in it to get the two pieces to transform consistently 

A\ = A {TT) a b + a~ 1 [LY} a b , 
A {TT) \ = 4>- Q A {TT) a b , [LY]\ = LY% , a = <t> 6 a , (43) 

so that the vector potential equation takes the form 

v b (a- 1 [LY\ b a ) = v b A b a (44) 

in terms of the barred variables, where the right hand side can be replaced using 
the supermomentum constraint Eq. (|41|l . This is the final improved conformal 
approach discovered by York |37j. The only obvious candidate for an extra 
initial data variable that transforms as a should is the metric volume factor 
g 1 / 2 or the lapse N, which has the same transformation properties when the 
Taub function is fixed. A convenient choice is a = 2N. This removes the 
previous ambiguity in the decomposition process (decomposing before or after 
the conformal transformation) and is much more satisfying. 

This additional factor of a inside the divergence operator corresponds to 
orthogonality of the splitting of the tracefree extrinsic curvature into transverse 
traceless and longitudinal parts under an inner product which includes the lapse 
as a multiplier of the spatial volume integration factor g 1 ^ 2 , i.e., in the combina- 
tion 

(4) g l/2 = Ng l/2 ( 

see the discussion after Eq. (4.6) of Pfeiffer and York |37j). 
This is easily seen by an integration by parts in which the extra factor of the 
lapse multiplying g 1 / 2 cancels the factor of the lapse multiplying the conformal 
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Lie derivative before that integration by parts, allowing the covariant derivative 
to flip over to the transverse traccless part yielding zero. 

The initial value problem takes place at one initial hypersurface, but if one 
has an evolving spacetime solution of the Einstein equations, one can see what 
this decomposition looks like in terms of the evolving spatial metric, which in 
the limit near the initial hypersurface takes us to the "thin sandwich data" 
approach |42| to the initial value problem using the spatial metric and its time 
derivative (g a b,9 a b) as variables instead of the spatial metric and the extrinsic 
curvature {g a b,K ab ), with 

dab = -2NK ab + £f3g ab , K ab = -^(g ab - £^g a b) , (45) 

giving the transformation of variables in both directions. One can then un- 
ambiguously translate the decomposition of the tracefree part of the extrinsic 
curvature into a decomposition of the tracefree part of the metric velocity, or 
the "conformal metric velocity" 

~gJ TF) = -2NA ab + [LPU , A ab = ^& t ( TF ' - [LPU) . (46) 
The proper time velocity is also useful 

Lb {TF) = N-% b W = -2A ab + N-^LPU . (47) 

The obvious next question to ask would be, how do we fix the two parts 
(transverse and longitudinal) of the tracefree part of the extrinsic curvature in 
either point of view? In the thin sandwich picture how is the potential vector 
field Y a related to the shift vector field (3 a l 

Given er and g ab , starting from any tracefree symmetric tensor C ab , one can 
remove its divergence to get a transverse traceless symmetric tensor 

A$ T) = C ab - a- 1 [LV] ab , V b {a~ 1 [LV] ab ) = V b C ab (48) 

which then determines the transverse traceless part of the barred extrinsic cur- 
vature by the conformal rescaling. The longitudinal part [LY] a& of the tracefree 
extrinsic curvature is then determined by the supermomentum constraint (|41() . 
from which the transverse term drops out, leading to 

V^Llf, = V b A b a = h/ a r + ] a . (49) 

The subtracted divergence part can then be combined with the vector potential 
term 

A ab = C ab + a- x [L{Y - V)] ab , (50) 
In the thin sandwich picture, we can make the identifications 
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Thus the shift is related to the original vector potential by 

j3 a = Y a - V a . (52) 

Suppose we take a spacetime sliced and threaded in zero shift gauge f3 a = 0. 
This implies that the vector potential Y a must equal the vector V a which gen- 
erates the divergence of the tracefree part of the conformal metric velocity. 
On the other hand one can choose the latter vector to be zero, insisting that 
the conformal metric velocity be transverse, which forces the shift to equal the 
vector potential. Both possibilities describe the two most useful spatial gauge 
choices for the Bianchi type IX spacetimes where the spatial diffeomorphism 
group is exactly the symmetry group SO(3, R) of the rigid body problem used 
as an analogy for understanding the spatial diffeomorphism gauge freedom of 
generic spacetimes by Fischer and Mardsen 02] • Zero-shift gauge corresponds 
to space-fixed coordinates in the rigid body problem, while the transverse gauge 
corresponds to body- fixed coordinates, with the spatial metric diagonalized ex- 
actly as for the moment of inertia tensor in the rigid body problem. 

One can also compare the vector potential equation (see Eqs. 144(1 . (|49[l ). 
which is equivalent to 

V b (-A h a + (t- 1 LY) = , (53) 
with the minimal distortion equation of Smarr and York |44j 

V b (-<rA b a + L(3) = . (54) 

These differ only by the overall multiplicative factor a inside the divergence 
operation. If a — 2N is spatially constant or if the Taub time gauge (constant 
a) is used so that the spatially covariant constant factor of N = ag 1 / 2 can pass 
outside the divergence operation, then these two equations coincide and the 
vector potential for the conformal metric velocity is also a minimal distortion 
shift vector field, which can then be used to change the spatial gauge to the 
minimal distortion gauge in order to minimize the time-rate of change of the 
magnitude of the conformal metric velocity 02] . This is exactly what happens in 
spatially homogeneous cosmology where the lapse is spatially constant |45l I46| . 

However, it may be appropriate |47| to alter the definition of the minimal 
distortion shift by including the lapse in the spatial volume integration in the 
variational principle in parallel with the new approach to the vector potential 
equation |37j (see their remark after their Eq. (4.6)), as well as using the proper 
time conformal velocity. This change in the variational principle from 

|^(TF)a b ,(TF) bQ -l/ 2w l Aw 2 Aw 3 (55) 

to 

J l^ a b °g^ b a Ng 1 ' 2 ^ A lo 2 A lo 3 (56) 
= J (-2A a b + 7V- 1 [L/3] a fc )(-2i fc a + N~ 1 [L(3] a b ) Ng 1 / 2 ^ 1 A u 2 A tu 3 
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leads to a net factor of TV -1 in the new variational integral, so that upon varia- 
tion with respect to f3 a and integrating by parts (provided the boundary terms 
can be ignored in a general spacetime), one has this extra factor of N~ 1 in the 
new minimal distortion shift equation, making it identical to the new vector po- 
tential equation for generic spacetimes, which is very satisfying, revealing a nice 
compatibility between the initial value problem decomposition and the dynami- 
cal evolution gauge-fixing that would be valid in general spacetimes. In the new 
minimal distortion gauge, the magnitude of the proper time derivative of the 
conformal metric is minimized over a thin sandwich between two nearby time 
coordinate hypersurfaces with the spacetime volume element as the measure 
over this 4-dimensional region. 

5 The Spatially Homogeneous Case 

The spatially homogeneous spacetimes ^5] have served a valuable role as a 
testing ground for general relativity. Restricting the variational approach to 
Einstein's equations to the minisuperspace of such spacetimes of a fixed Bianchi 
symmetry type |48| requires examining the boundary terms for troubles |49| . 
The same boundary term considerations can lead to a disconnect between the 
supermomenta as constraint equations and as generators of (symmetry compati- 
ble) spatial diffeomorphisms in this context |SU]. In order not to be distracted by 
the many cases which can arise, it suffices to consider the "trouble-free" Bianchi 
type VIII and IX spacetimes where the symmetry specialization does not intro- 
duce any such problems since spatial divergences of spatially homogeneous vec- 
tor fields automatically vanish, making all troublesome spatial boundary terms 
zero, while remaining general enough within the spatially homogeneous space- 
times to not suffer from other problems |46|. 

The conformal transformation is naturally identified with the usual decom- 
position of the metric into its unimodular piece and over all determinant. Using 
matrix notation 

g = e*?g, detg = l, g = . (57) 

Without introducing the barred/unbarred notation to complicate matters, we 
can associate the conformal factor with 

4 = e 2 ^ , = e*/ 2 , cf>~ 6 = e- 3& . (58) 

The unimodular part (conformal metric) is in turn decomposed into a diag- 
onal metric matrix and a linear transformation associated with the action of the 
natural (symmetry compatible) restricted diffeomorphism group on the spatial 
metric expressed in a computational frame composed of spatially homogeneous 
frame vectors 02] 

g = S T g (D) S , g (D) = e 2 ^ , det S - 1 , Tr^ = , = (3 ± e± , (59) 
where 9m\ = e 2 @ and (3 = /3 a l+f3 ± e± (sum implied over +, — ) decomposes the 
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logarithm of the the full diagonal metric matrix into its pure trace and tracefree 
parts using the Misner decomposition |161 I5T] . 

If we identify the space sections of the spatially homogeneous spacetime 
with the Lie symmetry group (simply transitive case), then the spatial diffeo- 
morphism freedom compatible with the use of invariant computational frames 
is the associated group of automorphisms and translations of that group, which 
induces the action of the linear automorphism group on the space of spatially 
homogeneous tensor fields. For the semisimple case of Bianchi types VIII and 
IX, there are only inner automorphisms, also called adjoint transformations, so 
the action of the linear adjoint group on spatially homogeneous tensor fields is 
relevant. This group is 3-dimensional (isomorphic to SO(2, 1) and SO(3) re- 
spectively) and when the spatial structure constants of the invariant frame are 
chosen to have diagonal form, it acts on the space of spatial metric matrices 
tranversally to the diagonal metric matrices and so can be used to decompose 
the metric matrix. S belongs to this group, which leaves invariant the spatial 
structure constants 

QC/~tl Q—lm Q—ln /-~rc fan\ 

>~>l l> mnd a*~> b — ab ■ V uu / 

The basis of the adjoint Lie algebra is 

[feal^c — ac i [fca 5 fc&] — C C a fcfc c , Sfc Q S 1 = kfrS^ a . (61) 

The metric parametrization g = S T g^S can be thought of as a transfor- 
mation from the initial gauge (usually taken to be zero-shift gauge) to diagonal 
gauge where the spatial metric is diagonal and the new spatial frame is 
£a' = EbS~ lb a - These new frame vectors must commute with the nonspatial 
computational frame vector ea> = d/dt — B a e a > in order to define a computa- 
tional frame, which leads to the condition 

SS" 1 =Cj a k a = B a 'k a -> B a ' = ib a (62) 

for the new shift vector field B a e a <, identifying its new components with the 
"autmorphism velocities." These enter the extrinsic curvature or metric veloci- 
ties in the following way 

-^g- 1 g = ^S- l [-p°i-h + ^ a K, a ]S , (63) 

where the matrix 

JC a = \{k a + e-^k T a e^) 1 (64) 

is just the symmetric part of the adjoint matrix k a with respect to the diagonal 
spatial metric. The second term in the metric velocity, containing the factor /3, 
is the transverse part of the conformal metric velocity, while the third term is 
its longitudinal part. 

If one starts with zero shift gauge, then the extrinsic curvature matrix is 
K = —(2N)~ 1 g~ 1 g, and one can interpret this choice of metric parametrization 
as pointing towards diagonal gauge (for which the conformal metric velocity is 
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transverse) , with the vector potential for the tracefree extrinsic curvature defin- 
ing the associated shift vector field for that new gauge. This was noted in the 
"old" conformal approach decomposition where the factor of a was omitted |46| 
leading to the outdated identification (3 a — 2N<fi~ 6 Y a instead of the much nicer 
result f3 a — Y a in the new improved approach. Thus the new approach to the 
supermomentum constraint and minimal distortion shift gauge finds complete 
resonance with the mathematical structure of spatially homogeneous dynam- 
ics, and suggests extending it to general spacetimes by redefining the minimal 
distortion shift condition. 

6 Conclusions and Acknowledgments 

Spatially homogeneous spacetimes have provided good feedback to the general 
theory in many ways over the past half century since they were born at the 
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more with Jimmy during my year with him. I had not yet learned the lessons of 
my long association with relativity in Italy that followed beginning the next year 
with my second postdoctoral position in Rome with Remo Ruffini. We are not 
in this game just because of the mathematical or physical interest it has for us, 
but also for the human relationships that it allows us to foster, especially those 
crossing national boundaries, a career perk that we have been very fortunate to 
have. 

Giorgio Ferrarese has helped keep alive this spirit of social ties among col- 
leagues that is typically Italian, especially with his tradition of Elba conferences. 
I was fortunate to be able to accompany Abe Taub on his last visit to Italy at 
the 1989 such conference (in menory of Carlo Cattaneo) where I recall a pleasant 
conversation with Yvonne on the beach. I am happy to have this opportunity 
again to see her in person and convey my appreciation for her long and fruitful 
career that has continued to cross paths with Jimmy's and which recently won 
them jointly the 2003 Heineman Prize for Mathematical Physics (announced in 
2002) awarded by the American Physical Society and the American Institute of 
Physics. In a way this helped bring Jimmy and me together again in Rio last 
summer by catalyzing his attendance at the Tenth Marcel Grossmann Meeting 
(2003) through my serendipitous intervention and subsequently reawakened my 
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